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ABSTRACT 

In this paper we suggest the existence in the central regions of spiral 
galaxies of collisionlcss, scale-free, rigidly rotating, sclf-gravitating discs with 
spiral symmetry. Such discs must be truncated at a finite radius , and they 
must be stabilized and rendered self-similar by a suitable halo. The halo and 
the rotating disc share the self-similar class and must form together to arrive 
at the suggested state. We make comparisons with the well-known rigidly ro- 
tating, Kalnajs discs; one of which is axi-symmetric and finite while the other 
is infinite and decomposed into spiral modes. We find the self-consistent, self- 
similar, distribution functions in one and two dimensions in a rigidly rotat- 
ing, collisionless system. In the case of two dimensions we deduce the self- 
consistency condition for discrete spiral arms. We give an estimate of the 
disturbance created in the halo by the presence of the disc, and argue that 
the halo itself should be close to self-similarity. A very weak cusp in the halo 
may be necessary. The necessary spatial coincidence of the halo results in a 
kind of disc-halo 'conspiracy'. Finally the disc equations are formulated in 
'spiral' coordinates, and the passage to an approximately discrete 'line spiral' 
is given as an example. Although in two dimensions the collisionless particles 
enter and leave the arms in non-linear epicycles, they move approximately 
parallel to the arms in the line spiral limit. The spiral pattern is however 
in rigid rotation. Aperiodic spiral arms are suggested wherein discontinuities 
may be coarse-grained to appear as collisionless shocks. 
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1 INTRODUCTION 



Axially symmetric, self-similar, thin, self-gravitating discs have been well studied; partic- 
ularly so in the works of Evans (1994), and of Evans and Read (1998a), which developed 
earlier work by Zang (1976). Such studies have revealed the surface-density/potential pairs 
and have deduced self-consistent surface distribution functions. Evans and Read (1998b) 
have investigated the stability of such discs by expanding the modes in logarithmic spirals. 

The conclusions regarding the stability of strictly scale-free discs were left somewhat 
uncertain until a summary paper by Goodman and Evans (1999) in which scale-free gaseous 
discs were examined. It was argued there that without boundaries, the normal modes of 
the scale-free disc were not properly defined. This was consistent with Evans and Read 
(1998b)who imposed a 'cutout' centre to study the normal modes. These modes tended to 
be neutral for bi-symmetric modes relative to the more rapidly growing one-armed or axially 
symmetric modes. 

The object of such time-dependent perturbation analyses is naturally to clarify the origin 
of the spiral structure of late type galaxies. Rather than attempt this discussion in this 
series of papers, we will instead bypass this question of origin by constructing non-linear, 
non-axially symmetric, stationary, uniformly rotating, spirals. In general non-linear stability 
requires a suitable stabilizing halo for equilibrium. The halo remains passive in our study 
and is thus artificial, but the existence of rigidly rotating structures in the centres of galaxies 
may nevertheless render its properties of interest. 

The stability of the rigidly rotating, entirely self-consistent Kalnajs disc was studied 
originally in Kalnajs (1972) and as part of a more general study by Vauterin and Dejonghe 
(1995). This latter work is very close to the present one, including the assumption that 
the disc is not isolated but is rather coupled to an indefinite passive halo, be it either of 
dark matter or of a bulge of stars or both. This halo is such as to make the disc potential 
harmonic regardless of the disc distribution function. The interesting result of this previous 
work for our purposes is that with the Kalnajs DF (as defined below) the uniformly rotating 
disc is found to be linearly stable in a sufficiently strong (relative to a maximum disc) 
constant density halo. No spiral structure develops. However they find this to be peculiar 
to the Kalnajs disc and NOT characteristic of uniformly rotating discs in general. In fact 
it is readily shown that the non-linear, non-axially symmetric distortions of a Kalnajs disc 
depend critically on both the central and edge cut-offs for equilibrium. 

Another work that uses surprisingly similar techniques is (Syer&Tremaine 1996). These 
authors also consider scale-free discs, but they replace the the collisionless matter by a 
barotropic, two-dimensional fluid. Their scaling constraints are necessarily somewhat differ- 
ent because of this, but nevertheless their /3 can be identified with our a — 5 (not to be 
indentified with our ad hoc use of (3 in section 3) and our £ is their 0. They also discuss 
(redone independently in this work!) the scaling of the Poisson integral as do we in the 
section on the Kalnajs disc. The self-similar case that they identify with — 1/4 is our class 
a = 5/4, while the special case /3 = is our case a = 1. Their weakly non axially-symmetric 
limit corresponds to our limit of small winding angle e/5, which we do not elaborate in this 
work. 

In the present work we treat both the non-linear spiral arms and the background axi- 
symmetric disc as self-similar (assuming the stabilizing halo) and collisionless. Thus we are 
taking the stars and possibly molecular clouds to be the dominant material component both 
of the disc and of the spiral arms. The resulting discs resemble locally infinite discs, but 
in fact they must be globally finite and stabilized as alluded to above and in more detail, 
below . Embedded in a real galaxy they are probably cut-off in a ring-like structure before 
the corotation resonance. 

We have no novel answer to the 'anti-spiral' objection (Lynden-Bell and Ostriker, 1967), 
other than to appeal to time-irreversible disturbances (e.g. 'fly-by' encounters) and dissi- 
pation. In one of our non-axi-symmetric solutions we shall in fact find that, discontinuities 



which imply shock structure of some type are required. It should also be noted that collec- 
tive effects in N-body systems (such as Landau damping) can lead to local perturbations 
which are soon amplified in their effects by Liapunov divergence. This leads to an effective 
irreversibility since virtually infinite accuracy in particle motion is soon required to time 
reverse the system. 

In this first paper we concentrate on the uniformly rotating disc, because it seems not 
to have been previously studied as a scale-free example (although the scale invariance of 
the Poisson integral was known; Kalnajs, 1971) and because it allows a straight forward 
introduction of the basic ideas. We shall extend the treatment subsequently to a more 
general scale-free disc. Such a constantly rotating disc does not correspond to the outer 
structure of any late type galaxy, although it is relevant to the inner structure that may 
display spirals and bars. It can also be applied to some irregular galaxies. 

Our treatment finds in most cases the surface-density/potential pair and the collisionless 
distribution function self-consistently by an 'ab initio' self-similar analysis. In section 2 
including various sub-sections we give the usual (e.g. Carter and Henriksen, 1991; Henriksen, 
1997, Le Delliou, Henriksen and MacMillan, 2011) self-similar formalism and deduce the 
general forms. Sections 3 and 4 deal with some axi-symmetric limits and make contact with 
the well known' Kalnajs disc' (Kalnajs, 1972; Binney and Tremaine, 1987, 2008). The non 
axi-symmetric Kalnajs discs whose potential follows from the decomposition of the Poisson 
integral into spiral modes, are also studied in these sections. Section 4 deals particularly with 
the halo structure required to permit stable, self-similar centres of finite discs. In section 
5 we construct various types of stationary, spiral arms, including discontinuous examples 
wherein shocks must appear. In section 6 a useful formulation of the problem in 'spiral 
coordinates' is given and then is used to study isolated spiral arms in a stabilizing halo. 
Finally we summarize and discuss our results. 



2 SELF-CONSISTENT, SELF-SIMILAR DISC EQUATIONS 

In this section we review the fundamental equations and deduce the self-similar forms for a 
rigidly rotating collisionless disc. We imagine it to be composed of stars and perhaps small 
molecular clouds. 



2.1 Basic Equations 

We work in cylindrical coordinates {r, <p, z} rather than in spherical coordinates because 
the algebra is simpler. However spherical coordinates are more appropriate if the disc is to 
be made compatible with external spherical mass distributions, and we use these in a brief 
discussion of the stabilizing halo below. The equations of the stationary thin disc are: 

Collisionless Boltzmann Equation in the Disc 

We use consistently a reference frame that is rotating with uniform angular velocity Q, so 
that both the azimuthal velocity v<j> and the azimuthal angle are measured in this reference 
frame. We designate the two dimensional distribution function (DF) by F(r,(j),v r ,v^) and 
the potential by $(r, (/)). We frequently use the notation d x for the operator d/d x where 
convenient. Other notations are standard. The equation is formally identical in spherical 
coordinates. 

v r d r F + —d<j,F + (^ + 2Q.V4 + Q 2 r - d r $)d Vr F - (^ + 2Qv r + -d+fyd^F = 0. (1) 



Laplace Equation Above the disc 



i<9 r (r<9 r $) + -^<9j$ + <9 Z 2 $ = 0. (2) 



Boundary Conditions for Laplace 

We designate the surface density of the disc by er(r, </>). We require for self-consistency 
the relations 

(<9$) z=0 = 2irGa, (3) 
and 

a = J F dvrdvtjy, (4) 

where F is a solution of equation (1) using the potential $(r, 0, 0) = $ (r, 0). 

So far we have considered the disc as infinite in radius, but such a self-similar, uniformly 
rotating disc can not be in equilibrium 1 . Thus there is a radial boundary that must be 
considered, plus a surrounding halo that is required for equilibrium. We will discuss these 
matters in detail below, but it suffices here to say that for a finite disc a halo can be contrived 
that permits the central regions of the disc to be self-similar as described here. The required 
halo does not differ greatly from a uniform density spherical mass distribution, which is also 
of similarity class a = 0. 

This completes the set of basic equations. 



2.2 Self-Similar Analysis 

We transform the basic equations to 'scaling invariant variables' according to the usual 
analysis (see e.g (Le Delliou, Henriksen & MacMillan 2011); (Henriksen 1997); (Henriksen 
& Widrow 1995);(Carter & Henriksen 1991)). The transformation takes the form: 



Independent Variables 

Here 5 is the (reciprocal) spatial scaling factor and the mass scaling factor is related to it 
as /i = 35, because G is invariant and the time scaling factor a = 0. The latter is imposed by 
the invariance of Q under the scaling. We retain the remaining scaling parameter 8 explicitly 
below because we can sometimes simplify results in the coarse or fine grained limits as 5 — > oo 
or 5 — > respectively. Apart from this possible use, we may set 5 = l/r where r D is an 
arbitrary fiducial radius. Thus the scale-invariant variables R, Z, £ are found from (Carter 
& Henriksen 1991) 

Sr = e SR , ^ = (j) + eR = (f)+(e/5)\n6r, 6 z z = Ze SR . (5) 

We observe that the non-axisymmetric self-similar variable is a logarithmic spiral in vari- 
ables {r, <ft} with rd(f)/dr = —e/8 as the tangent of the 'winding angle' equal to ninety degrees 
minus the pitch angle. The extra scaling parameter 5 Z can be useful in strict boundary-layer 
problems where it might be taken much smaller than 5 in order to ignore the z dependence 
in the potential, or much larger to accentuate it. We will normally take it to equal 5 below. 



Dependent Variables 



1 I am indebted to a referee for emphasizing this point 



We transform the independent variables to functions of the scale invariants {R, £, Z} ac- 
cording to the dimensions of each quantity while taking R as the Lie scaling group parameter 
(see e.g Henriksen and Widrow (1975)) to write 

$ = e 2SR V(R,£,Z), F = e~ SR P(R, £, Y) 

a = e SR X(R,£), v = e SR Y. (6) 

These follow our standard procedure where for example F = exp (2a — 45 + /x) in terms of 
time, space and mass dimensions and scales. However a = for fixed Q and \x = 35 — 2a to 
conserve Newton's constant under scaling, and so finally we obtain the specified F = e~ SR P. 

It will be verified below in the section on Kalnajs discs that the Poisson integral for the 
potential in the plane may be wholly expressed in terms of these variables. 

We proceed by a direct substitution of these dependent and independent variables into 
our set of basic equations. We obtain the equations that the scale independent quantities 
must satisfy in the forms: 



Self-Similar Boltzmann 

Our substitutions yield 

Y R d R P + (5Y lf> + eY R )d i P+(5(Y^-Y^ + 2QY (t) + Q 2 /5-ed^-d R ^-25^)d YR P 
-(25Y^Y R + 2{lY R + 5d^)d Y ^P = 5Y R P. (7) 

This equation is easily solved by characteristics (the self-similar particle trajectories), which 
will be carried out below. 

Self-Similar Laplace 



We find that our substitutions yield for this equation 

+■ (^) 2 ) e$v + i e 



Sftf + — + + 4^ = 0. (8) 



Self-Similar Boundary Conditions 

Using 5 Z = 5 the boundary conditions become 
2ttGS = 5(d z V) , (9) 
and 

r P dY R dY$. (10) 



It is important to realize that we impose the self-similarity ultimately by requiring the 
scale-invariant functions P,\& and S to be independent of the group parameter R. This has 
already been enforced explicitly in the self-similar Laplace equation. 

Is the Self-Similarity Physical? 

A referee has properly questioned the physicality of an infinite self-similar disc on the 
grounds that the outer regions of such a disc always dominate the local gravitational field. 
This would lead locally to an infinite gravitational force directed outwards, in stark contra- 
diction to the assumed equilibrium. 

This correct argument implies that the self-similarity can only apply to a truncated 
disc of radius 'R' at sufficiently small r, where one might not expect the disc radius to be 
physically able to destroy the self-similarity. In fact one finds that for most surface density 



profiles a suitable halo is still necessary to stabilize the disc and to impose the self-similar 
behaviour. A non axially symmetric disc is somewhat more suited to the self-similarity, but 
in general the halo and the finite disc must be self-similar together. We shall discuss these 
points in detail below. 

In the next sub-section we study the scale-free Boltzmann equation in order to see how 
this symmetry restricts the normal Jeans' theorem. 



2.3 The Disc Scale-Free Distribution Function 

In this section we study the solution to equation (7), with and without axial-symmetry. 
Proceeding generally, it is easily found that the characteristics of equation (7) yield together 

§ = - 2SY « £ ' (12) 

where ds is along the characteristic and £ is the scale-free energy given by 

These latter two equations may be combined to yield the general value for the scale-free 
distribution function in the form 

p = Wr>> (14) 

where K is strictly constant. It may depend on a characteristic constant if one can be 
found, but no such constant is available in the strictly self-similar equations. 

We may translate this result into non-scaled variables by using the transformations given 
in equations (6). In this way we find 

F = Wr^ (15) 

where the characteristic energy in the rotating frame is E = v 2 /2+Q— Q 2 r 2 /2. It is important 
to note that if we write explcitly the upper limit in energy as E Q , then in strict self-similarity 
this must be of the form kQ 2 r 2 where k is a real number. This is due to a constant energy 
being incompatible with the a = self- similarity class. 

It is instructive to relax the self-similarity by allowing P{R) and ^(-R). This does not 
change the two equations (11), (12) above (although it does change the Poisson equation) 
but one now obtains in addition that Yr = dR/ds. This means that these latter two equations 
may be integrated directly to give £ = Ee~ 25R and P = Fe SR . We have in these expressions 
identified the two constants of integration to be the the characteristic energy E, and the 
distribution function F, also constant on a characteristic. These identifications follow from 
the scalings of equations (5) and (6). 

Consequently by substituting these integrals, equation (14) becomes 



lEn-El 1 / 2 ' 



(16) 



Unlike the situation in the self-similar limit, one can take the energy E a to be constant. 

Basically in this non-self- similar case, this says only that F = F(E) as by the Jeans' 
theorem. However it does suggest a special case where K(E) is constant. This proves to be 
useful in our discussion of the Kalnajs disc. 

Should the solution we seek be axi-symmetric, the azimuthal velocity on the characteristic 
satisfies from equation (7) 



dY 

^± = -2Y R (5Y (P + n), (17) 

which is not integrable in the self-similar case. However by dropping the self- similarity and 
using again Yr = dR/ds we obtain the integral 

5Y^ + Q = 6 2 3z e~ 25R , (18) 

where once again we have used equations (6, 5) to identify the constant of integration as 
essentially the specific angular momentum 5 2 j z . This implies a two integral DF which we 
may take in the suggestive form 

n K(E, E/(Qj z )) 



(19) 



By letting E/(Qj z ) be universally constant we return to equation (16), after which we may 
impose the self-similar form by taking K(E) constant and setting E a = kQ 2 r 2 . 

We recall that the energy in the rotating frame E is related to the energy in the inertial 
frame Ej by E = Ej — ttj z , so that taking E/(Qj z ) constant ensures that the ratio Ej/Qj z 
is also constant. This latter constant is equal to unity in WKB density wave theory (see 
e.g. equation (6.86) of Binney and Tremaine 2008) and in fact also for the multi-poles of 
electromagnetic radiation (Jackson, 1999). The definitions are slightly different in each case, 
but a ratio of unity seems to reflect an essential relation between angular momentum 'flux' 
and energy 'flux' in a wave field. A uniformly rotating disc may be the 'near zone' of such 
a wave field. 



3 THE 'KALNAJS' DISCS 

In this section we first compare the rigidly rotating Kalnajs disc to the (finite) self-similar 
disc. We note that although the surface density at small radius is not self-similar, both the 
potential and the self-consistent DF are. Subsequently we discuss the Kalnajs (Kalnajs 1971) 
disc/potential decomposition into logarithmic spirals and show in what sense this approach 
is self-similar. These are not rotating discs in strict self-similarity since a ^ 0. In fact we 
find in subsequent work (noted below) that the corresponding DF is isotropic in an inertial 
frame. 



3.1 Rigidly Rotating Finite Disc 

For the usual Kalnajs disc it is convenient to drop the scale- free formulation. Formally with a 
fixed angular velocity and a fixed radius both the temporal scaling a and the spatial scaling 
S are zero. The only remaining self-similar 'motion' is through the dependence on £, but this 
does not concern us in this section since we do not include spirality. 

This disc is defined by taking the flat limit of an axially symmetric MacLaurin spheroid 
of radius a to find the surface density 

2 

a = a /(l- r -). (20) 
v or 

For this reason it is also referred to as a 'MacLaurin' disc. The remarkable property of this 
disc is that in the plane, and for r < a, the potential has the same form as does the infinite 
self-similar disc namely 

Q 2 r 2 

$ = \ , (21) 
where the circular angular velocity is given by 



2a 

This disc (Kalnajs, 1972) has been discussed recently in (Binney& Tremaine 1987) and 
(Binney &Tremaine 2008) and an explicit form for the potential everywhere has also been 
given recently by (Schulz, 2009). In the latter work the potential found in the plane for r < a 
can be easily verified by using the Bessel function version of the density/potential pairs. 
The effective potential in the rotating frame is 

*J - ^^r 2 , (23) 

It is easy to verify that for r < a a consistent solution requires Q c > Q so that $j > 0. 
We note that because of symmety in v<p in the DF (16), the average value of this quantity is 
zero. This means that the average circular velocity in the inertial frame is = Qr, giving 
Q a physical interpretation. 

We can now close the loop by evaluating the surface density from the integral (4) using 
the DF (16) with general E . This procedure yields with $j > and E Q > 



I / Q2 _ Q2 

a = AnK\J E Q - $, = 4irK^(E - r 2 ), ( 24 ) 

and this must equal cr ^Jl — r 2 /a 2 for consistency. Insisting on the equality requires 

n 2 - n 2 

E = (25) 
which is always greater than <3>j as is necessary; and the further condition on K as 

** K = 1 J° n2 • (26) 



This demonstration of self-consistency was given already in (Kalnajs 1972), but it is some- 
times discussed only in the inertial frame (e.g. (Binney & Tremaine 1987), (Binney & Tremaine 
2008)). We have presented the argument here because the self-similar approach in a rotating 
frame gives a simpler, more inevitable description. For example, the radial and azimuthal 
dispersions are obviously equal in this frame and the mean rotational velocity is immedi- 
ate. We also emphasize the connection with the distribution function of the scale-free disc, 
and in fact with the distribution function for radial orbits in a logarithmic potential (Le 
Delliou, Henriksen & MacMillan 2011), which is also of the form (16) in 3D (the Fridmann- 
Polyachenko (Fridmann & Polyachenko 1984) type). 

We emphasize once again that to obtain the self-similar disc from the DF (16) one may 
take E D = k 2 (Q 2 — Q 2 )r 2 /2 in equation (24). There then follows the scale-free surface density 
as 



'n 2 - n 2 



a ss = AnKVW^lsj r. (27) 

But as we shall see below, the self-similar disc must also be truncated at finite radius, which 
causes the potential to differ from the quadratic form. Even worse, such an isolated disc can 
not be in gravitational equilibrium. This must be corrected back to the self-similar form by 
an appropriate halo mass distribution for the self-similarity to be physically realizeable. The 
necessary halo is itself of some physical interest however, since rigidly rotating structures do 
exist in the bulges of spiral galaxies. 

It is interesting to demonstrate as we remarked previously that assuming Ej/(Qj z ) is 
constant suffices to pick the Fridmann-Polyachenko (or Kalnajs) DF from the family of axi- 
symmetric possibilities (19). For when the DF is a function only of E > we have the 
functional form (Kalnajs, 1976) for the DF as 



F ^ = -hwy (28) 

Recalling equations (23), (25), (26) and (20) this expression yields 

F <*'> = Je^ww" (29) 

and hence agrees with the functional form F(E) of (16). The argument can not be used in 
this form for the self-similar limit of this DF, because of the dependence of the upper limit 
on r and hence on $j. 



3.2 Self-Similar Spiral Decomposition 

We start with the Poisson integral over the surface density for the disc potential, and apply 
a self-similar analysis. We do not now assume uniform rotation so that scaling in time a is 
generally not equal to zero. The usual dimensional analysis and Lie group formalism yields 
(note that —25 + fi <— 5 — 2a because of the invariance of G under the scaling) 

$ = ^(£) e 2 ( 5 - Q ^, a = S(Oe {5 ~ 2a)R 

Sr = e SR , v = e {s - a)R Y. (30) 

We substitute these into the Poisson integral to obtain a relation between the scaled 
(also known as 'reduced') potential and surface density as 

*(£) = — / du ^ / / (31) 
o J-oo Jo ^(2 cosh(n) - 2 cos(</>' - </>)) 

where 

u = 5(R'-R), ^' = (f)' + eR' 

= 1 -2a, a = j. (32) 
2 o 

The case studied by Kalnajs (ibid) and presented in (Binney &Tremaine 2008) is in fact 
the only strictly self-similar case (see also (Syer&Tremaine 1996)), since then — and 
the integral converges with the infinite limits. Otherwise finite limits are necessary which 
introduce a dependence on R. 

We have nevertheless introduced this slight generalization because to the extent that R 
is small compared to say the radius of a finite disc A, the upper limit is then approximately 
independent of R. Then for convergence > — 1/2 (a < 3/2) is required. This allows another 
class of approximately self-similar discs to be found, but without the z dependence of the 
potential we are unable to impose self-consistency with the DF (see below). An infinite disc 
with a central 'cut-out' would require < 1/2 (a > 1) while a finite disc with a central 
cut-out allows general 0. The integral is difficult in such cases, however. 

The Kalnajs case requires = and so we infer the corresponding self-similar 'class' 
(Carter & Henriksen 1991) as 

•-H- (33) 

Using these values in equation (30) reveals that \P and £ are identical to the so-called 
'reduced quantities' used in (Kalnajs 1971) and (Binney &Tremaine 2008). Kalnajs (Kalnajs 
1971), referring to earlier work by Snow, may have been aware of the Lie self-similar nature 
of the transformation to logarithmic radius and reduced variables. However the idea was 
used neither generally nor explicitly as we have done here. The scaling was discussed in a 



similar but more ad hoc fashion in (Syer&Tremaine 1996), and their parameter f3 = 1/4 is 
equivalent to the self-similar limit above. Our a = 1 class would correspond to their class 0, 
which we see would give a divergent potential unless the disc is finite. We do not discuss this 
case here as it is not uniformly rotating, but later in this series we study it as the general 
type of 'Mestel disc'. 

To proceed we seek, following (Kalnajs 1971), to express equation (31) in terms of spiral 
modes. We write this following (Binney &Tremaine 2008) 2 ) as 

ny] poo r 2w p im(C+^u) 

UJ-T 1 ^/ du dC, . * (34) 
d J -°° Jo J2cosh(u) - 2cos(() 



where £(£') = E m e !m ^ and £ = <j>' — <$■ The double integral, say N(m, e/5), may be evaluated 
according to (Kalnajs 1971) as 



2 



^^)-^Tfrr#F- (35) 

We have verified this result numerically to less than one half percent in particular case with 
m = 2 and e/5 = 1 . 

Although this result is an elegant derivation of a potential/density spiral pair in the 
plane, it does not help us to produce self-consistent spirals since we do not have the z 
dependence. This can be resolved by finding the solution of the Laplace equation above the 
disc. We leave this treatment mostly to a subsequent work however, since such a disc does 
not fall into the category of uniformly rotating discs. 

We have carried out an analysis similar to that given above for the uniformly rotating 
self-similar disc, which will be reported in detail elsewhere. It is perhaps worth recording 
here nevertheless that we find the DF that corresponds to the disc of self-similar class given 
in equation (33) to be 

F = K(E) 2 . (36) 

This gives another example where the upper limit in an integration over energy must be 
compatible with the self-similar class. This implies that an upper limit E a must have the 
form kr~ x l 2 . One readily finds that the integral of the DF (36) over velocity space with such 
an upper limit does give a oc r _3//2 as required. 

In the next section we treat the axially symmetric scale-free rigidly rotating disc. It is 
an approximation to a finite disc given an appropriate stabilizing halo in a sense that we 
study in some detail. We may, when the disc is not wholly composed of spiral arms, add this 
azimuthally uniform background to the non-axi-symmetric case. The non-axially symmetric 
potential thus acts as a non-linear perturbation on the axially symmetric background. 



4 AXIALLY SYMMETRIC, SCALE FREE, RIGIDLY ROTATING DISC 

In this section we discuss the problem of stabilizing a self-similar disc by a passive halo of 
stars (bulge) or dark matter. For this to be possible the disc must be finite. This requirement 
is not removed even by taking the 'disc' to be composed solely of a small number of spiral 
arms. 

We start by calculating the potential at radius r in an axi-symmetric disc of finite radius 
R that has the expected surface mass density a = E 5r (S constant). We do this by using 
the result found in equation (2.155) of (Binney &Tremaine 2008) to write 



2 There is an unimportant difference since to maintain formal self-similarity £ must be the independent variable so that m also 
multiplies the term in radius 



\/r 2 — a 2 J a \f- 



" a — a 2 



The inner integral is straight forward and after some rearrangement we obtain 
r l R dx ( o, / x 



$ d = 2GZ5R 2 / ; [x 2 In 7— ~ V^x 2 , (38) 

where x = a/R. 

To obtain the gravitational acceleration we proceed similarly by using equation (2.157) 
as given in (Binney &Tremaine 2008). This allows us to write the gravitational acceleration 
in the disc at radius r as 

iGSSR2 <" R dx , 1 L H *) + , : ) . (39) 



dr r Jo /((r/RV - t 2 \ \ 1 -i- , /C\ - ^ 



((r/R) 2 -x 2 ) V l + y/(l-x 

The integrals in these two expressions are plotted in figure (1) as a function of r/R. 

We observe from these figures that the gravitational force is outwards for r/R < 0.65 and 
inwards nearer the edge of the disc. Thus the disc is unstable to forming a very thin ring at 
the transition radius! Moreover the potential shows no tendency to behave self-similarly even 
though we have forced it to be zero at the origin. There is a small region between 0.65 and 
about 0.85 where the potential can be fitted with a quadratic law, but this is nevertheless 
in the context of a globally unstable disc. 

In order to stabilize this behaviour while maintaining self-similarity we must require a 
halo to provide a potential $/i(r) at the disc such that 

$ & (r) + $ d (r) = V a (6r) 2 . (40) 

The disc will in principle ensure that the total halo potential is not spherical 3 but assuming 
that the halo dominates, it should remain roughly spherical. Applying the Laplacian to 
equation (40) yields roughly the halo density ph as 

4,G Ph = «>». - 2 -'^-% (41) 
r dr dr A 

The last two terms in this expression give the spherical approximation to the non-uniform 
core formed from the undisturbed halo. We may refer to them as i Deltarhd h and their form 
in units of r = r/R is shown on the lower left of figure (1). Note that for technical reasons 
r is written in the discrete form 0.01 + 0.05j where j runs from 1 to 19. We see that the 
presence of the disc establishes a very weak cusp in the otherwise uniform (cored) halo. 

This evaluation of the disturbed halo density is a crude approximation, provided that 
the disc contribution to the total halo potential taken spherical remains small. This will not 
be the case as the centre of the system is approached arbitrarily closely, because of the weak 
cusp. However the cusp is weaker than logarithmic (which is compatible with the variation 
found in numerical simulations and some theoretical estimates (Henriksen 2007)). Given that 
uniformly rotating barred and spiral structures do exist in the centres of galaxies, it may 
be worth taking semi-seriously this form of the cusped halo. When the non-axi-symmetric 
structure is inside a ring, even the outwardly increasing surface density may be relevant. 

In order to estimate the potential disturbance produced in the halo by a truncated self- 
similar disc near its centre and not too far above the disc, we take the self-similar forms 
literally. Then with axial symmetry and E/Qj z constant, the DF is the self-similar result of 
equation (14). The Laplace equation (8) becomes simply 

d 2 z m + 4^ = 0, (42) 



3 It is difficult to know this disc potential for z ^ but we find a near-disc approximation below. 
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Figure 1. The figure at upper left shows the truncated self-similar disc potential &d('>~) in units of 2GT,SR 2 as a function of 
r = r/R. We sec that its gradient changes from negative to positive at r = 0.65 so that the disc is unstable to ring formation. 
This is seen more directly on the plot of the potential gradient d^^/dr at upper right. At lower left is shown the required 
perturbation to a uniform cored halo, assuming that the halo remains spherical. The zero of the potential has been placed at 
the origin as is the case for the self-similar potential. 



whence 

* = Acos2Z + Bsin2Z, 

where A and B are constants. In the disc 
A-Q 2 /(25 2 ). It seems that this should be > 
mean circular velocity in the inertial frame, 
for the axi-symmetric surface density 

X a = 47rK(y/E o -yJVj{0)), 



(43) 

we have that the effective potential ^j(O) = 
since otherwise gravity could not balance the 
In this case the boundary condition (10) yields 

(44) 



that is explicitly 



/ — / fi 2 

S a = 4ttK(\J S q - \J(A- — ). (45) 

The energy £ Q is simply an upper bound to the self-similar particle energy. Its physical 
correspondant scales with r 2 . 

The remaining boundary condition (9) takes the explicit form 

B = ^S a . (46) 
o 

For fixed E a and £ Q the constants A and K appear to be degenerate. However if, as 
for the Kalnajs disc, we require gravity to give the required centripetal acceleration for the 
mean circular velocity, then A = Q 2 /(25 2 ) and \I/j(o) = 0. The potential is then determined 
everywhere as 

= — cos2Z + ^-£ a sin2Z, (47) 
2o z o 

or in terms of non-scale invariant quantities 

^ n 2 r 2 2z „ . 2 . 2z 

$ = cos h 7rGS a 5r 2 sin — . (48) 

2 r r 

When the disc is stabilized by a halo potential A is equal to (5r) 2 ty a where \l/ a is due to 
the halo plus disc. For completeness we note that the self-consistent DF amplitude is given 
in general by equation (45). 

The potential (48) is pathological as r — > or as z — > oo but it is well-behaved on radii 
with z/r constant, except for the divergence at radial infinity. Similar problems arise in 
spherical symmetry, including radial divergence and a logarithmic singularity at any point 
on the polar axis. It is likely that, when using this infinite solution to approximate a finite 
disc of radius R, we will require z <C R and so the pathology will not appear off the axis. 

Equation (48) does not allow consistently for the finiteness of the disc. It merely assumes 
finite behaviour in the self-similar form. Nevertheless its form suggests that there is only a 
weak dependence on the spherical angle tt/2 — 9 = arctan(z/r) near the disc. 

A uniform density halo obeys the a = self-similarity 4 including the quadratic potential, 
so a uniformly rotating disc and a uniform halo are of the same self-similar class. The halo 
is not expected to be uniformly rotating so that it must acquire its self- similarity from the 
boundary condition at the disc. However we have seen that the stable existence of the disc 
requires a distortion of the halo that is not self-similar. Provided that this distortion is 
small, the common self-similarity class plus the need for disc stabilization provides a kind 
of 'disc/halo conspiracy' on the scale of the galactic nucleus. 

From this section we conclude that the self-similar forms do have some application to the 
central parts of finite discs. In the next section we proceed therefore to give a description of 
the spiral structure of a uniformly rotating self-similar disc. 



5 SPIRAL/BAR STRUCTURE IN THE RIGIDLY ROTATING DISC 

All that is new in this section stems from the solution for the non axi-symmetric potential 
that is found from equation (8). This solution we can expect to be superimposed on the 
axisymmetric solution of the last section. In effect we are treating a non-linear example of 
the Lin-Shu hypothesis (see e.g. (Binney &Tremaine 2008)) where the pattern speed is the 
uniform value Q and the spiral form is stationary in the rotating frame. 



4 In 3D the density dimension co-vectord p = (0, —3, 1) — > (0, 0) so that the self-similar density is uniform. 



Equation (8) is easily separable and readily yields a solution in the real domain. Unfor- 
tunately such a solution is not strictly periodic in the plane. It is a Laplace transform mode 
in z rather than a Fourier mode in the plane. This leads to discontinuities in the surface 
density and in the potential. Nevertheless, because the analysis is transparent relative to the 
periodic discussion that follows, and because the discontinuities may well imitate line shocks 
in the plane, we discuss this case in the next sub-section. Such shock structure is probably 
required to select preferentially either a trailing or a leading spiral, and so to escape the 
'anti-spiral' theorem. 

These shocks would be in collisionless matter and although they are usually associ- 
ated with plasmas they have been suggested for the rotating gravitational N-body problem 
since the collective interactions are similar (Fridmann & Polyachenko 1984),(Fridmann et al 
1996),(Binney &Tremaine 2008),(Mouhot& Villani 2009). 

The last reference in this series argues that non-linear Landau damping leads to small 
scale filamentation in phase space, which might signal the transition to a new equilibrium 
state. The first and second references in this list deal specifically with rotating systems and 
suggest the transition from soliton to shock by non-linear dissipation. The thickness of the 
shock is in any case finite (less than an epicyclic radius in the soliton mechanism and perhaps 
slightly greater than some effective Jeans length). 

Moreover there is an associated gravitational potential gradient that leads to a change 
in the potential across the dissipative region. Seen macroscopically after coarse graining, 
this implies a jump in the potential across the shock discontinuity. This process may be 
a reasonable description of the collisionless transition to a spiral arm, and it could either 
initiate or be initiated by a hydrodynamic shock. We employ this shock hypothesis in the 
next section, but subsequently we return to periodic arms without shocks or distinction 
between leading and trailing examples. 



5.1 Aperiodic Spiral Structure 

We separate the equation (8) according to \l/ = S(£)C(Z),and we choose ( oc exp (— k\Z\) 
(where k > 0). Moreover to remain as close to a pure planar mode as the equation permits 
we can set for integer m 

k 2 ^k 2 m = m 2 (l + ^)--^ Tj (49) 

1 + 

as this choice makes S(£) periodic except for the damping (equation(50) below). This im- 
position requires e 2 /5 2 > 1 for m — 1. Thus the one-armed spiral 'shock' (see below) is 
relatively tightly wound, as this requires the 'winding angle' arctan |e/<5| to be greater than 
45°. For m = 2 or greater the only requirement on the winding angle is that e 2 /5 2 > 0. 
The solution for each 'mode' of the scaled potential becomes now 

V m (Z, = A m exp [ - (k m \Z\ + 2e/ / 0] cos (m£ + 6 m ), (50) 

1 + e z /o z 

where A m and m are real constants. This solution has the form of a damped harmonic 
oscillator and so is path dependent with regard to the variation in £. This is a kind of 
irreversibility disguised as non-integrability. 

To complete the solution we must satisfy the boundary conditions (9) and (10) at the 
disc. The first condition gives the form of the surface density as 

£(0,£) = E„ + E m = E o -MV27rG)* m (0,0 = 

E a - 7^ k m A m exp [ - ( - ^Jj 5 2 0] cos ( m £ + Qm), (51) 



where the potential is evaluated at Z = 0. The surface density E a describes an axi-symmetric 
component. This surface density corresponds to the net stabilized axi-symmetric potential 
(Sr) 2 V a . 

The second boundary condition gives an expression for the surface density of collisionless 
matter. We assume that the net axi-symmetric potential is positive in the inertial frame for 
stability, and require that \& m (0,£) + \P a = Q 2 /(25 2 ) + * m (0, £) > 0. In the rotating frame 
we suppose the net potential energy \l/ a + \l/ m (0,£) — Q 2 /(25 2 ) = ^ m (0, £) of the spiral 
disturbance to be negative in the peak density regions. 

We might very well assume that the peak energy is positive. With the negative assump- 
tion the perturbed gravitational force is directed outwards in the rotating frame and so 
particles might be expected to cross the arm on the outward leg of their 'epi-cyclic' motion. 
With the positive assumption they would cross the arm on the inward leg. This suggests 
that the arm has a forward facing shock while the negative assumption implies a backward 
facing shock. We proceed with the negative assumption, but it should become clear how to 
deal with the positive energy assumption. 

We obtain from equations (10) and (14) that in regions where the spiral energy is negative 

f a Y, a + / m S m (0,O = 47rA^ m (0,O| 1/2 , (52) 

where /j < 1 is the fraction of the disc material {i} that is collisionless. The origin of 
(1 — must be in another form of matter, almost certainly gas and dust in a spiral 

galaxy. It is of course possible that f a = f m - 

We substitute for the spiral surface density on the left of this last equation using the 
appropriate part of equation (51). In this way we obtain a condition for self-consistency for 
each value of m in the negative energy region. Strictly one should only apply this condition 
when it contains the disturbed potential summed over all m, but in view of the nearly 
pure modal behaviour of at least the 'grand design' spiral galaxies we explore each mode 
independently. 

In regions of the disc where the spiral energy is positive the right hand side of equation 



(52) must be replaced by 47rA(-y/£^— y \l/ m (0, £)) as P er our previous discussion. An essential 
difference in this non-axially symmetric disturbance is that £ D can be a function of £, without 
destroying the self-similarity. In fact in order to pass smoothly through the zeros of \& m (0, £) 
it should have the form k 2 cos (m£ + m ) in the positive potential regions. 

The self-consistency condition in the negative energy regions becomes explicitly 



fn 



-k m A m exp [ - ( ) f J cos (mf + 9 r 



2ttG y l + e 2 /S 2 ' 



+ fa^a 



= 4 7 rA|A n exp[-( £ ] cos (m£ + 9 m )| 1/2 . (53) 

1 + e z /o z 

Dimensionally we see that everything is coherent so long as A m has the dimensions of squared 
velocity. 

There are various ways to satisfy this equation that we now consider. It is apparent that 
one choice of modal constants {A m , m } will not allow this condition to be satisfied for all 
£, for constant K, S a and fa. It is possible to see this condition as an equation for the single 
fraction /(£) = fa provided that its value remains less than one. However this leaves much 
arbitrariness in the problem. In order to obtain discrete spiral arms we shall require that 
fi < 1 be known values on the arm 5 . This will determine the arm. Since then the condition 
(53) can only be satisfied at discrete values of £ by fixing the modal constants. The two 
dimensional nature of the analysis therefore implies that particles in the background pass 
to and from the discrete spiral arms. 



5 Any fixed values would serve, but ultimately they are likely to be fixed by observation. 



One does not need to satisfy the second boundary condition outside of the spiral arms 
but it is likely to apply to the collisionless part of the disc. One must then treat the equation 
(53) as an equation for /(£) once we add the information that f a — fm — f ■ This means 
that / will be discontinuous across the arms. 

To illustrate the possible solutions, we consider first the one-armed mode with m = 1. 
Let us place a discrete spiral arm at £ = by choosing A\ > and ©i = it, so that we have 
a positive spiral density on the background according to equation (51). We recall that in 
this case we are restricted to e/8 > 1, so that the one-armed winding angle is always greater 
than 45°. Then equation (53) determines A\ from 

A\ + (2 V 2 - v 2 jf)A x - r? 4 = 0, (54) 
where 

8n 2 GK 2 2vrGS a . . 

v " "m"' 77 " "SfcT' (55) 

and v, t] have the dimensions of velocity. This equation always has one positive root when 
v 2 > ^rff 2 . All of the modal constants are now fixed by choosing / < 1. 

Starting at the positive or leading 'side' of the £ = spiral and incrementing <fi by 2ir at a 
fixed R brings us to the negative or trailing side with a density lower by a factor exp (— 2nb). 
We have defined for convenience 

6 "TT^W (56) 

which is always less than or equal to 1. There is thus a backward facing 'shock' type discon- 
tinuity from the trailing to the leading edge of the spiral. The spiral potential has the more 
negative value by the factor exp (2nb) after a 2ir increase in at fixed R. This is really a 
coarse-grained description of the potential gradient in a collisionless layer or shock. 

There is another type of solution that gives a set of discrete spirals. The same value for 
Ai given by equation (54) will apply on discrete spirals £ s that satisfy 

e" fe6 cos£ s = l. (57) 

This equation has an infinite number of negative roots which, after a least negative root, 
oscillate in £ s between close to —270° and —90°. We have converted the radian measure of £ s 
to degrees (modulo 360) in order to illustrate the separation in angle at a fixed radius. For 
example if we choose e/S = 2 so that b = 0.8, than the first root (after zero) is —66°. 843, and 
the succeeding two roots are —271°. 297 and —89°. 893 (modulo 360) respectively. These 
spirals are shown in figure (2). The fourth root spiral is too small to be seen on this scale 
that is chosen so that r = l/5 = r o at0 = Oon the £ = spiral. 

In this way we have constructed a set of one-armed spirals, each of which will have a 
density jump from the leading to the trailing side of magnitude for trailing spirals. 

The reverse is true for leading spirals. 

One can also use the continuity of the collisionless fractions to determine them in terms 
of A m by writing equation (53) at £ = and at £ = 2n and solving for the /j. This yields 
formally 




(58) 

These must be less than one but A m remains arbitrary otherwise. 

For the two-armed spiral with m = 2, the same array of solutions is possible as were 
found for the one-armed spiral above. We can once again start with an arm at £ = and 



Figure 2. The figure shows on the left a polar plot of one-armed aperiodic spirals that are solutions to equation (57) with 
e/<5 = 2. In descending order of size the spirals are at £ = {0, —1.1666, —4.7350, —7.8521} respectively. One sees that they are 
separated by about tt along the y = axis. The two-armed spiral on the right has discrete arms at f = (right spiral at 
y = 0) and at § = tt (left spiral at y = 0). The lower figure shows the same two-armed spiral when e/5 = 0.2, giving a bar-like 
structure. In all cases the stellar motions should be along the arms for physical consistency. 

satisfy the equation (53) by insisting that A 2 satisfy the equation (54) wherein k\ replaced 
by &2. For symmetry with a second arm at £ = tt, we require a jump in the density between 
the leading and trailing sides of magnitude e bn . Once again the reverse is true for leading 
arms. 

The main difference between the two-armed spiral and the one-armed spiral is that we 
are allowed to construct bars in this case by allowing e/5 — > 0. On the right-hand side of 
figure(2) we show a two-armed open spiral of this type having e/5 = 2 together at lower left 
with a nearly bar-like structure having e/5 = 0.2. 



Just as in the one-armed example we can construct a system of spiral arms located where 



e~ Ks cos2£ s = 1, 



(59) 



although the constant on the right could be any constant c that would then be absorbed 
into A 2 . The effect is to delete the arm at £ = 0. 

The roots of this last equation behave similarly to those found in the one-arm case. 
Using once again the value b = 2, there is a first root after the zero root at £i ~ —139°. 12 
and then the roots oscillate (modulo 360°) between approximately —45° and —225°. In fact 
they converge on these latter two values. Since each value of £ s found in this manner has its 
counterpart at £ s — n, these convergent pairs are the natural two-armed spiral. To ensure 
this we must remove the extraneous roots at zero and the first root £i. This can be done by 
setting c > 8 as the constant on theright of equation (59). In this way we arrive naturally 
at a two-armed spiral structure. Each arm will have a density jump of e a ^ from the leading 
to the trailing side (and the reverse for leading arms). 

One might think that this 'bipolarity' that arises in the one-armed case and in the two- 
armed case is an argument for the ubiquity of two-armed spirals. However, similar behaviour 
occurs for m = 3 where the interval between the convergent roots is ~ 120° as would be 
expected. The principal difference as m increases is that the strength of the 'shock' is reduced 
(to e 26?r/3 for m = 3) for the same winding angle. Thus, although the bi-polarity is not 
favoured by the system, the definition and hence the visibility of the arms may be maximal 
in the two-armed mode (after the true one-armed system). In this connection we note again 
that b(e/8) is a maximum (of one) at e/8 = 1 so that the strength of the discontinuity 
is reduced as the winding increases towards tight spirals or decreases towards bars. The 
interesting feature of this solution is that, because the oscillating roots become rather dense 
as the sequence continues, it is possible to imagine a true 2D PDF that defines the arms in 
this type of solution. 

In fact a true 2D solution to the self-consistency condition (53) also exists even with / 
constant, if the amplitude of the DF is K(£). This is degenerate with a variation in /(£). 
But such a variation requires £ to be a characteristic constant of the CBE, which in turn 
by equation (7) requires SY^ + eY R = 0. This is the condition that the particles in the disc 
in the rotating frame move purely backwards and forwards along the spiral direction, since 
Y^/Yr = rdcj)/dr = —e/8. Unhappily this requirement reveals Y R and to be dependent 
variables, so that the Boltzmann equation (7) is no longer correct. One of Yr, Y^, must be 
integrated out subject to the constant £ condition. This procedure is much more transparent 
in spiral based coordinates that we introduce below. 

The constants f m and f a could be determined also by continuity as we did for the one- 
armed case. 

Having studied spirals with shock type dissipation, we turn in the next section to consider 
the purely periodic, modal solutions to equation (8). Although these have a familiar modal 
structure, they do not offer a consistent means of avoiding the anti-spiral theorem. 



5.2 Periodic Spiral Structure 

In this section we seek periodic solutions of equation (8) by writing 



from which we must extract the real part. One must sum over all positive and negative integer 
values of m in order to describe an arbitrary periodicity, but we shall content ourselves 
with looking at pure modes. These modes, we recall, are in the rotating frame and are 
superimposed on the axi-symmetric background. By substituting the assumed form (60) 
into equation (8), we obtain an equation in the complex plane for Q{Z) as 



* = e- iim °((Z), 



(60) 



+ [4(l-^l)-(l + l!) m 2 ]C = O. (61) 

The potential for negative m and leading spirals is thus the same as that for positive m and 
trailing spirals. 

This harmonic equation is readily solved formally, but it can be tedious to extract the real 
part. It is much easier to consider individual modes. The solution in cylindrical coordinates 
has singularities at large z and small r, which is an apparent argument against a modal 
analysis in these coordinates. However we recall that the potential is everywhere well-behaved 
on a cone generator of constant Z, which is effectively the same behaviour as is found 
in spherical-polars. Indeed an analysis in spherical polars does not change the essential 
behaviour of the stellar disc. The CBE is identical in the plane to the one we have used even 
when written in rotating spherical polars. In any case we know that this analysis holds only 
with a stabilizing halo and a truncated disc not too far from the plane. 

Starting with the one-armed structure when m = 1 we obtain as a solution for ((Z) 6 



Ci(Z) = d sin (Vl - 4ie - e 2 Z) + C 2 cos (VI - 4*e - e 2 Z). (62) 

The constants C±, C 2 are complex in general, but since in their complex form they only add 
a phase to the sine and cosine functions, we take them to be real for simplicity. The real 
part of Z) still possesses a formidable form, but in the plane it becomes simply 

ite(*i(0,0) = C 2 cos£. (63) 
The derivative d^i(^, Z)/dZ is slightly more involved at the plane as 



Re = f U/2 v/(e 2 + 9)(e 2 + 1) + 6 - 2e 2 (cos£) 



e 2 + 9)(e 2 + l)-6 + 2e 2 (sinf) ), (64) 
which is related to the total surface density through 

^ = * + *g**{w\- (65) 

We assume once again that the spiral disturbance has negative energy in the rotating 
frame and that \l/ a = Q 2 /(25 2 ) so that equation (52) applies. Then the self-consistency 
equation becomes on using equations (63) and (64) 



fx [C 1 (y2 v /( e 2 + 9)( e 2 + l) + 6-2e 2 cose - ^(e 2 + 9)(e 2 + 1) - 6 + 2e 2 sin^)] + fJ^^ 
16vr 2 KG 
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C 2 cos£| 1/2 . (66) 



Even with a complex C 2 that would add a phase constant to £ under the cosine, there is 
no solution for the constants that would apply for all £ with K and fi constant. We proceed 
therefore as in the previous section by setting the spiral arm at £ = with fi = f a = f 
fixed. Otherwise we suppose that equation (66) gives /(£) in the negative energy region(and 
the corresponding equation in the positive energy region). 

For the arm at £ = 0, the constant C± > follows by setting the desired surface density in 
the arm from equation (65) given the winding angle. Subsequently, C 2 follows from equation 
(66) given K and /. There is no restriction of winding angle and no need for shock-type 
discontinuities. Both trailing and leading arms are possible. Once again collisionless particles 
must enter and leave the arm. 



In the subsequent equations of this section t <— e/6 for simplicity, unless otherwise noted 



For the two-armed mode with m = 2 the solution for ^2^, Z) is formally (the constants 
are now peculiar to this mode) 

^ 2 = (C x e 2 ^V&^ z + c 2 e- 2 ^V^ z ^j , (67) 

from which the real part at Z = follows as (we take C\ and C 2 real) 

ite(¥ 2 (0)) = (C t + C 2 ) cos2£. (68) 

Similarly some calculation shows that 

Re(— ) = (Ci - C 2 ) (6+(e) cos2£ + 6_(e) sin2£) . (69) 
We have set for brevity 

b ± = [2e v / ?T4±2e 2 ] 1/2 . (70) 

For large e, that is for very tightly wound spirals, 6_ — > 0. 

The self-consistency condition can now be written for negative energy spirals as 

MCa - C7 a )(6+(e) cos2£ + 6_(e) sin2£)] + / a ^^ = ^^V|(Ci + C 2 )co8 2£|. (71) 

Once again we see that the arms are discrete. They can be taken to lie at £ = and at 
£ = 7r, while again f2 = f a = fon the arms. Otherwise the equation determines /(£) for the 
inter-arm regions as above. 

There are two constants C\ — C 2 and C\ + C 2 to be determined. The self-consistency 
condition (71) at £ = or n will determine one for fixed K and /. However if the underlying 
axi-symmetric disc is given, so that S a is known then C\ — C 2 may be determined by giving 
the total arm density. 

It is worth noting that in the rotating frame there is no net flux of particles because 
of the isotropic nature of the Kalnajs DF. Thus particles enter and leave the spiral arms 
following what would be essentially epicyclic orbits in the inertial frame in the linear theory. 
In the multi-arm case particles may travel between the arms. 

We turn in the next section to an analysis based on the spiral structures as coordinate 
lines. This enables us to easily pass to a treatment where only motion along discrete spiral 
arms is present. 



6 SPIRAL ARM BASED ANALYSIS 

We begin this section by transforming from plane polar coordinates {r, 0} to spiral coordi- 
nates {£, £}, where £ is measured as length along the spiral and £ plays the role of orthogonal 
coordinate. This requires a careful representation of the Boltzmann equation although the 
Poisson equation changes very slightly. 

6.1 Equations in Spiral Coordinates 

The unit vectors locally parallel and orthogonal to a spiral arm defined by 

£ = + ^ln5r, (72) 

may be taken respectively as 



and 



h - ^ (74) 

These lead us to introduce the corresponding velocity components 

vi = e e -v, v $ = e 5 ■ u, (75) 

and these define in turn the length along the spiral d£ = Vi(£/)dt and an orthogonal velocity 
r£ = r<fi + (e/S)r = v $ + (e/S)v r . The first definition integrates to 




(76) 

while the second expression can be written using the definition of as 

" e= (i+W (77) 

Note that we also have by definition that vg = £■ 

We are now in a position to follow the usual route to the two dimensional CBE in a 
uniformly rotating frame by expressing the disc particle Lagrangian in terms of {£, £, £, £} 
in this frame, hence passing to the Hamiltonian, and finally to dF/dt = dF/dt + [F, H] = 
([F,H] is the Poisson bracket). 

The resulting canonical form does not hold the velocities {ve,v^} constant in the partial 
derivatives since the canonical momenta {p£,p^} are the independent variables. Hence ap- 
propriate transformations of the partial derivatives (including to (OF / d£)p^) must be made 
to obtain the following form 

d t F v e d e F + V^v&F + ( tl±2^i + v l _ d #) g p 




(78) 

This equation is general in a uniformly rotating frame. One sees that it reduces to equation 
(1) as e/5 — > (when the time dependence is removed for stationary rotating structures) as 
it should. The presence of the winding angle as a parameter opens the way to an expansion 
of the CBE either 'near' (in parameter space) bar-like structures (e/5 — > 0) or 'near' ring-like 
structures for which e/5 — > oo, but we do not pursue this here. 

The Poisson equation above the disc can be written simply in {£, £, s} coordinates, where 



'1 + e 2 /5 2 z, as 

jdt{id&) + $ + <9> = 0. (79) 

This equation together with the previous equation may be used as usual for a two dimensional 
disc to impose self-consistency. The results do not differ of course from those found above, 
although the description is perhaps more natural. However whenever a single spiral arm is 
to be discussed, as is the case in this section, an infinite number of 'modes' of this equation 
(to form a delta function in £) are required to satisfy correctly the boundary condition at 
the arm. We shall see that it is better to use the Poisson integral in such a case since then 
we simply bypass the difficult solution for the potential above the plane in favour of that in 
the disc. 



6.2 Discrete Spiral Arms 



We wish in this section to consider the case where in the rotating frame the particles move 
only parallel to the spiral arm. This is true one dimensional motion in our rotating spiral 
coordinates. Realistically, the motion should extend over a finite region in the £ coordinate. 
We find that although such finite structures are not so violently unstable as is the finite 
two-dimensional self-similar disc, they do still require weak adjustments to the halo in order 
to ensure self-similarity and stability. 

We assume that the collisionless distribution function on a spiral line can be written as 

F-FMv ^Vf-U , (80) 
where Sd{x) is the Dirac delta function. An integration of equation (78) over and £ yields 

vid t F t + - d e $) d Ve F e = 0, (81) 



1 + 



where Fg and $ are evaluated on the spiral £ = £ s . 

To give the spiral arm a modestly finite extent in £ we can interpret the delta function 
according to 

s (t t \ j ■ y ^ e(e-6)-e(e-6A 

but replace the strict limit by a finite interval [£1,^2]- Here Q(x) denotes the Heaviside step 
function, equal to 1 when the argument is positive and zero otherwise. 

We are now in a position to carry out a self-similar analysis for the DF on discrete spirals 
in a uniformly rotating frame of reference. For line masses the uniformly rotating (a = 0) 
self-similar mass per unit length A oc r 2 since d\ = (0, —1,1) — > (0, 2). We write, following 
the usual technique with a = 0, 

S£ = e 5L 5s = Ze 5L 

F e = P e e SL $ = e 2SL <H v t = Y e e 6L (83) 

and substitute into equation (81) to obtain 

Y L d L P e - (5Y? + d L m + 25V - y )d YL Pe = -5P t . (84) 

We recall that only when = do we have strict self-similarity, but even allowing this 
not to be the case the following conclusion holds. The characteristics of the latter equation 
allow us to conclude that (ds here is an increment along the characteristic) 

^ = -25Y L £, (85) 

where £ = Y?/2 + * - tt 2 /25 2 , and 

= -trtPt- (86) 

ds 

Together these two equations imply that the one dimensional DF on a spiral is 
P-£ = K\£\ 1/2 , (87) 
or equivalently in physical variables 

F = K\E\ 1/2 . (88) 

Only when there is no self-similarity does Yl = dL/ds, which allows equation (85) to be 
integrated directly to give £ exp 2L = constant = E = v 2 /2 + ^ — Q 2 £ 2 /2. In that case, just 



as we found in our discussion of the two dimensional DF at the beginning of this paper, we 
might have K = K(E) and so revert to the Jeans' theorem. 

Equation (88) reveals a satisfying progression in the self-similar DF with dimension for a 

system in uniform rotation. In one dimension it is oc \J\E\ as found here, in two dimensions 

it is oc 1/ yJ\E\ as found earlier in this paper, and in three dimensions it is oc 1 / (^J \E\) 3 . 
This latter result may be inferred from (Le Delliou, Henriksen & MacMillan 2011) for the 
isotropic limit when a = 0, although the energy would have to be written in the uniformly 
rotating frame as above. The basic result as a function of a in the maximally coarse-grained 
steady limit may be found in earlier papers such as (Henriksen 2007) and references therein. 

For strict radial infall in spherical symmetry, the self-similar DF is also found to be oc J\E\ 
(Henriksen & Widrow 1995), (Le Delliou, Henriksen & MacMillan 2011) 

We can imagine this collisionless spiral arm to be superimposed on the usual axi- 
symmetric, truncated self-similar disc, for which both the density and the potential (quadratic 
after stabilized by the halo) are known. This may be made in part of collisionless matter as 
above. However here we do not allow a constant interchange of particles between the arm 
and the disc 

Just as in the two dimensional case, one should verify to what extent a truncated self- 
similar discrete arm can be stabilized by an appropriate halo. The potential, given the 
forms (83), satisfies the same equation above the plane as (8). The boundary condition on a 
discrete spiral is difficult to apply with the modal solutions of this equation. We use instead 
the general self-similar Poisson integral (31). 

We write the surface density in the form (recall the more physical expression of the delta 
function) 

a = MLzifl A ( r ) = £(£)(<f r ), (89) 

where the linear mass density has the self-similar form A = A(<5r) 2 . Together these last 
expressions give 

£(0 = (A5)M£-6)- (90) 
For more than one arm there would be a sum over the various £ s (each with its own K s , 
giving an effective K(£)) but for present purposes we take only one arm at £ s = 0. It should 
be clear how to generalize to two or more arms from the following. 
The physical potential now follows for a = from equation (31) as 

= _^ e -a/,*/' (9i) 

v2 J-oo J cosh (x' — x) — cos (£ s — — (e/8)x') 

We have taken a finite disc of radius R ma x = 1/5 and x = 5R = In Sr. The physical potential 
is written rather than the scaled potential since the potential, being long range, might 
scale differently from the local quantities such as the surface density and the DF. In such 
a case the self-similar Poisson equation would have to start again from equation (79), but 
since we do not use this equation we can ignore this requirement. 

We have taken the physical delta function to evaluate the integral over <p' at the average 
value —{e/S)x' + (£i + G)/2. The average (£ = (£i + ^)/2 = £ s ) f° r simplicity for small 
£ 2 — £2- However this creates a logarithmic singularity due to a rigorous delta function that 
must be avoided. This can be done formally by keeping £ s and £ distinct. In our present 
example we set £ s = and discuss only one arm. For several arms, A = A(£ s ) and there 
would be a sum of terms over £ s , which label is manifestly an integral of the motion for the 
CBE. 

To find the potential on a spiral £ = constant we constrain in equation (91) to equal 
£ — (e/S)x This gives finally 



-0.3 



Figure 3. The left hand figure shows the run of the potential in r at small r for different angles. Labelling the curves 
at r = +0.1 and starting from the top the angles are {7ir/&, 3n/2, 7r/3600, 7r/4, 57r/4, tt/2, it, 3tt/4} . The right hand figure 
shows the run of the potential in £ at different radii. Labelling the curves at £ = 5, the radii in units of the disc radius are 
{0.1,0.05,0.001,0.005,0.01}. The constant value of the potential at r = 0.001, £ = tt/3600 equal to 5.226684385 has been 
removed in all cases so as to show better the variations. 



GA a/2), f 1 e^^'dx' 

$ = -=e- [l/2)x / =. (92) 

v2 -'-oo ./cosh (x' — x) — cos ((e/5)(x — x') — £) 

We observe that there is a logarithmic singularity on the line £ = at x' = x for all x since 
the difference of functions under the square root goes to zero as x' — x) 2 . However this may 
be avoided by retaining £ 7^ in the argument of the cosine. The integral is readily evaluated 
numerically for £ as small clS cL milliradian. 

We have evaluated the integral (92) numerically both as a function of r (recall that 
x = \ri8r) and of £. We see the results in figure (3). 

We conclude that the single armed disc is weakly unstable (outward force where the 
gradient in r is negative), which has to be corrected by a suitable variation in the halo. 
More interestingly we see that for r < 0.1 (in units of the disc radius) $ is rather constant 
to within a few per cent (recall that a constant has been removed from the potential of 
the order 5.227). The scaled potential \P is certainly not independent of r and must be 
corrected to the quadratic form. In deriving the Boltzmann equation we have assumed that 
a halo maintains the self-similar harmonic potential. There would only be a weak correction 
to a uniform density halo given the small variability at small radius, but the necessary 
'conspiracy' remains. 

Physically, we have found in this section a way of constructing discrete rigidly rotating 
spiral arms that are stabilized by a halo. The self-consistent DF is plausibly given by equation 
(88) although we can not study the self-consistency condition without knowing the potential 
above the plane. The particles move purely along these arms in the rotating frame. There is 
no net flux of particles along the arms because of the isotropy of the DF, so a particle going 
out at one stage is returning along the arm after a period in the rotating frame. 



7 DISCUSSION AND CONCLUSIONS 

In this study we have looked at spiral structure in finite thin discs from the perspective of 
stationary self-similarity in a rigidly rotating frame of reference. It is necessary to appeal 



to a suitably modified uniform density halo (dark or otherwise) to ensure the existence and 
the stability of such a structure. 

We have examined the nature of the disequilibrium in a finite disc that has the linear 
surface density profile. In the case of axi-symmetry we have estimated the perturbation to 
the uniform density halo required for stable self-similarity. A kind of disc-halo conspiracy 
is required in general. In later sections we find the form of the potential in the plane of a 
discrete spiral arm, which also requires corresponding perturbations to the halo. 

We found that the 'Kalnajs DF' of equations (14) and (15) is the same as the self-similar, 
uniformly rotating, disc DF for collisionless matter in the rotating reference frame. This is 
so even though the surface density is not self-similar, although the potential that it inher- 
its from a MacLaurin spheroid is. We were able to fit this DF into a dimensional sequence 
wherein it characterizes a two dimensional distribution, while F oc l^ 1 / 2 is the correponding 
distribution in one dimension and F oc |i?|~ 3 / 2 is the equivalent distribution in three dimen- 
sions. When one relaxes self- similarity one obtains the slightly modified Jeans' distribution 
functions in equations (16) and (19). We remarked that the DF (19) is nevertheless reduced 
to (16) if one supposes that E/Qj z is constant, as is the case for spiral waves. 

The equivalent truncated self-similar disc has the same DF and harmonic potential as 
does the Kalnajs disc, but it has a different upper limit in energy space. This limit is 
determined by the self-similarity and is quadratic in radius, being proportional to Q 2 r 2 . 
This is the only difference with the Kalnajs disc, given the existence of a halo that renders 
the disc potential harmonic. The self-similar disc can not exist in isolation. 

In this same section we reviewed the scaling symmetry of the Poisson integral for a disc 
potential and so found the general scaling for finite discs. In the particular case due to 
Kalnajs (Kalnajs 1971), wherein the potential is expanded in spiral modes, we saw that the 
similarity class in the sense of (Carter & Henriksen 1991) was a = 5/4. 

In section 4 we studied the corrections to a truncated self-similar disc potential that 
are necessary to obtain a harmonic potential. This is a measure of the disc-halo conspiracy 
required. Continuing with the self-similar disc including a surface density linear in radius, 
we derived briefly the potential structure above the plane of the disc. Its main function is 
to provide a background medium for the non-linear spiral arms discussed in the following 
sections. 

In section 5 we studied first aperiodic, non-linear, spiral modes. Each of these require 
spiral discontinuities that should correspond physically to collisionless 'shocks', perhaps 
including the possibility of hydraulic jumps (Martos&Cox 1998) . This built-in mechanism 
for dissipation probably justifies focussing on trailing arms, although formally both cases 
are possible. The self-consistency condition (53) requires discrete spiral arms for constant K 
and constant fractions of collisionless matter. There is a mode in which the arms may come 
in clusters. 

It is interesting to note in this connection that the one-armed structures are constrained 
to be more open than a winding angle of 45°. Since the surface density increases outwards, 
the one-armed structures might resemble more one-sided bars than rings. There is no such 
restriction for two-armed structures so that they might well be mistaken for rings at radii 
corresponding to the physical end of the disc. 

We turned next in this section to study true periodic spiral perturbations and looked 
at the individual modes. In this case there are no discontinuities, but other considerations 
remain much the same as in the aperiodic case. The self-consistency conditions (66) and (71) 
lead once again to discrete arms on using the self-similar DF with K and fractions f m and 
f a constant. By using the self-consistency condition for a single mode we have assumed that 
each mode dominates in turn. We avoid that assumption finally by constructing isolated 
spiral arms that necessitate only an additional stabilizing halo. 

To assist this latter construction we changed coordinates from rotating plane polars to 
rotating 'spiral coordinates'. The Boltzmann equation in these coordinates has its uses for 



discs, but we employed it here to describe a model for a spiral arm in which the motion is 
wholly parallel to the arm in the rotating frame. The motion is such in the inertial frame 
that a particle will execute a loop in which both its outward and inward motions are parallel 
to the arm. By using an approximate delta function for the surface sensity, the arm may be 
constructed to have a finite extent in the spiral coordinate. It is thus a bundle of spirals all 
with the same winding angle. 

Strictly one must avoid concentrating all of the mass on a line since then the arm potential 
diverges and can not be corrected by a finite halo. Accepting this, the potential is found to 
be only weakly varying at small radius by evaluating the scaling Poisson integral. Thus the 
halo may readily impose the harmonic potential without too much distortion. In this case 
the distribution function on the arm is given by equation (88). 

Our treatment has not permitted us to say anything about time dependence, except that 
such structures would not exist in the absence of a suitable halo. The stability is most readily 
achieved if the halo also is close to the self-similar class of the disc in three dimensions. This 
is a kind of 'disc-halo conspiracy'. Physically these structures can only be taken literally in 
the central regions of finite discs embedded in halos. However this may not be unreasonable 
in the central bulge regions of spiral galaxies. The pattern speed of the spirals is one of rigid 
rotation, but the collisionless particles generally move relative to the arms. This motion is 
either in and out of the arms in two dimensions or along the arm in the narrow arm limit. 

Finally one may well ask why should such an apparently contrived system be of interest 
in reality? It does appear that uniformly rotating discs exist in the bulges of spiral galaxies, 
but they might of course be of the Kalnajs type. The only distinguishing feature would be 
the linearly increasing surface density of the self-similar disc. Moreover it is clear that the 
halo and the uniformly rotating disc would have had to form in a coordinated fashion. 

The origin of the kind of conspiracy required may well be the coincidence of both the 
uniform density core and the uniformly rotating disc with the a = class of self-similarity. 
The formation of a dominant core provides the potential necessary for the equilibrium of the 
self-similar disc. Material falling onto the plane from a growing core would tend to produce 
a linearly increasing surface density provided that it also obeys the self-similar symmetry. 
That is, since there can be no fixed scale, the mass must fall from a height that increases 
linearly with r with a density that is everywhere constant. Moreover the thickness of the 
slice Ar associated with the mass at r, although small, must also increase linearly as kr for 
small k. Thus as we have often concluded, the self-similar disc must be truncated at a finite 
radius R such that kR «r. 
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